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at an invariant mass of W = 1125 MeV and four-momentum trans- 
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a Rosenbluth separation of the transverse and longitudinal cross sec- 
tions was performed. An effective Lagrangian model was used to ex- 
tract the 'axial mass' from experimental data. We find a value of Ma = 
(1.077 ± 0.039) GeV which is (0.051 ± 0.044) GeV larger than the axial 
mass known from neutrino scattering experiments. This is consistent with 
recent calculations in chiral perturbation theory. 
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1 Introduction 



There are basically two methods to determine the weak axial form factor 
of the nucleon. One set of experimental data comes from measurements of 
(quasi) elastic (anti) neutrino scattering on protons [TJ, deuterons and other 
nuclei ( Al, Fe) J3]|| or composite targets like freon and propane [§|§ . In 
the (quasi)elastic picture of (anti) neutrino- nucleus scattering, the uN — * /iN 
weak transition amplitude can be expressed in terms of the nucleon electro- 
magnetic form factors F\ and F 2 and the axial form factor Ga- The axial form 
factor is then extracted by fitting the Q 2 -dependence of the (anti)neutrino- 
nucleon cross section, 

= A(Q 2 ) T B{Q 2 ) {s-u) + C{Q 2 ) (s - uf , (1) 



in which Ga{Q 2 ) is contained in the bilinear forms A(Q 2 ), B(Q 2 ) and C(Q 2 ) of 
the relevant form factors and is assumed to be the only unknown quantity. It 
can be parameterised in terms of an 'axial mass' Ma as Ga(Q 2 ) = G\a(0)/(1 + 
Q 2 /Ml) 2 . 
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Fig. 1. Axial mass Ma as extracted from (quasi)elastic neutrino and antineutrino 
scattering experiments. The weighted average is Ma = (1.026 ± 0.017) GeV, or 
(1.026 ± 0.021) GeV using the scaled-error averaging recommended by Ref. |l0[ . 

Fig. [I] shows the available values for Ma obtained from these studies. Refer- 
ences reported severe uncertainties in either knowledge of the incident 
neutrino flux or reliability of the theoretical input needed to subtract the back- 
ground from genuine elastic events (both of which gradually improved in sub- 
sequent experiments). The values derived fall well outside the most probable 
range of values known today and exhibit very large statistical and systemat- 
ical errors. Following the data selection criteria of the Particle Data Group 
fllDj , they were excluded from this compilation. 



2 



Another body of data comes from charged pion electroproduction on protons 



TT|— pL7|| slightly above the pion production threshold. As opposed to neutrino 
scattering, which is described by the Cabibbo-mixed V — A theory, the extrac- 
tion of the axial form factor from electroproduction requires a more involved 
theoretical picture. 
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Fig. 2. Axial mass Ma as extracted from charged pion electroproduction experi- 
ments. The weighted average (excluding our result) is Ma = (1.068 ± 0.015) GeV, 
or (1.068 ±0.017) GeV using the scaled-error averaging [|n]]. Including our extracted 
value, the weighted scaled-error average becomes Ma = (1.069 ± 0.016) GeV. Note 
that our value contains both the statistical and systematical uncertainty; for other 
values the systematical errors were not explicitly given. SP: soft-pion limit, DR: anal- 
ysis using approach of Ref. pTl, FPV: Ref. pi BNR: Ref. pj. 



The basic result about low energy photoproduction of massless charged pions 
can be traced back to the Kroll-Ruderman theorem ||18|| , extended to virtual 
photons by Nambu, Lurie and Shrauner |19j, who obtained the 0{Q 2 ) result 
for the isospin (— ) (see Ref. [pSj] , p. 29 for notation) electric dipole amplitude 
at threshold 



0,Q 2 



1 
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0(Q 3 ) ,(2) 



where k v is the nucleon isovector anomalous magnetic moment, = Ga(0) 
is the axial coupling constant, and /„- is the pion decay constant. 

In the following years, improved models were proposed Hl6| , p0| -p2]] , most of 
them including corrections due to the finite pion mass n = m^/M. The values 
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of the axial mass were determined, within the framework of the respective 
model, from the slopes of the angle-integrated differential electroproduction 
cross sections at threshold, 



i) 



6 dG A (Q 2 



Ga(0) dQ 2 



Q 2 =0 



12 

Ml 



(3) 



The results of various measurements and theoretical approaches are shown in 
Fig. ^. Note again that references [0,0 were omitted from the fit for lack of 
reasonable compatibility with the other results. 

Although the results of these investigations deviate from each other by more 
than their claimed accuracy, the weighted averages from neutrino scattering 
and electroproduction give quite precise values of the axial mass. Comparing 
the average values of the two methods, one observes a significant difference 
of AM A = (0.042 ± 0.023) GeV, or (0.042 ± 0.027) GeV using the scaled-error 
averaging. 

Chiral perturbation theory (xPT) has recently shown a remarkable and model- 
independent result that already at 0(Q 2 ), the NLS result of eq. (0) is strongly 



modified due to pion loop contributions f25|| . These contributions effectively 
reduce the mean-square axial radius, 

(rl) (rl) + ^ (l - P) • (4) 



The loop correction in eq. (^) has a value of — 0.046 fm 2 , which is a —10% 
correction to a typical (r A ) = 0.45 fm 2 . Correspondingly, the axial mass Ma — 
\/T2/ (r 2 ^) 1 / 2 would appear to be about 5% larger in electroproduction than 
in neutrino scattering, in agreement with the observed AMa. The aim of 
the present investigation was to determine Ma from new, high precision pion 
electroproduction data and thereby help verify whether this discrepancy was 
genuine. Since both the energy and momentum transfers were too high to allow 
for a safe extraction of £'04- > these data were analysed in the framework of an 
effective Lagrangian model with the electromagnetic nucleon form factors, the 
electric pion form factor and the axial nucleon form factor at the appropriate 
vertices [0,0. 



2 Kinematics of the experiment 



The differential cross sections for ir + electroproduction on protons were mea- 
sured at an invariant mass of W — 1125 MeV and at four- momentum transfers 
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of the virtual photon Q 2 = 0.117, 0.195 and 0.273 (GeV/c) 2 . For each value 
of Q 2 , we measured the scattered electron and the outgoing pion in parallel 
kinematics at three different polarisations of the virtual photon, enabling us 
to separate the transverse and the longitudinal part of the cross section by 
the Rosenbluth technique. Table [I] shows the experimental settings. 

Table 1 

Experimental settings for the p(e, e'-7r + )n experiment. The angles are measured with 
respect to the electron beam axis. 

Setting (e) Q 2 E c E' e 6 C 6 n p n 

[GeV 2 /c 2 ] [MeV] [MeV] [°] [°] [MeV/c] 



0.834 


0.117 


855 


11 


587.35 


-27.93 


39.31 


188.84 


0.500 




510 


11 


242.35 


-58.22 


28.31 




0.219 




405 


11 


137.36 


92.96 


-18.41 




0.742 


0.195 


855 


11 


545.79 


-37.72 


38.27 


209.62 


0.437 




585 


11 


275.89 


66.67 


-28.03 




0.229 




495 


11 


185.79 


93.45 


-20.12 




0.648 


0.273 


855 


11 


504.55 


46.83 


-35.82 


228.00 


0.457 




690 


11 


339.55 


65.27 


-29.37 




0.259 




585 


11 


234.55 


89.60 


-21.90 





3 Experimental setup 



The measurements were performed at the Institut fur Kernphysik at the Uni- 



versity of Mainz, using the continuous-wave electron microtron MAMI |2"6 



The energies of the incoming electron beam ranged from 405 to 855 MeV, and 
the beam energy spread did not exceed 0.16 MeV. The 15 to 35 /xA electron 
beam was scattered on a liquid hydrogen target cell (5x1x1 cm 3 with 10 
Havar walls in settings with the photon polarisation parameter e = 0.437, 
0.648, 0.457 and 0.259, and on a 2 cm-diameter cylindrical target cell with 
50 /iin Havar walls in all other settings) attached to a high power target cool- 
ing system. Forced circulation and a beam wobbling system were used to avoid 
density fluctuations of the liquid hydrogen. With this system, luminosities of 
up to 3.2 ■ lO^cm-V 1 (32MHz//ib) were attained. 
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The scattered electrons and the produced pions were detected in coincidence 
by the high resolution (5p/p ~ 10~ 4 ) magnetic spectrometers A (SpecA) and 
B (SpecB) of the Al Collaboration (27|. In settings with e = 0.834, 0.500 and 



0.742, electrons were detected with SpecB and pions with SpecA, and vice 
versa in all other settings. The momentum acceptance Ap/p was 20% and 
15 % in SpecA and SpecB, respectively. Heavy-metal collimators (21 msr in 
SpecA, 5.1 or 5.6 msr in SpecB) were used to minimise angular acceptance 
uncertainties due to the relatively large target cells. 

A trigger detector system consisting of two planes of segmented plastic scin- 
tillators and a threshold Cerenkov detector were used in each spectrometer. 
The coincidence time resolution, taking into account the different times of 
flight of particles for different trajectory lengths through the spectrometer, 
was between 1.0 and 2.6 ns FWHM. 

In each spectrometer, four vertical drift chambers were used for particle track- 
ing, measurement of momenta and target vertex reconstruction. Back-tracing 
the particle trajectories from the drift chambers through the magnetic sys- 
tems, an angular resolution (all FWHM) better than ±5 mrad (dispersive and 
non-dispersive angles) and spatial (vertex) resolution better than ±5 mm (non- 
dispersive direction, SpecA) and ±lmm (SpecB) were achieved on the target. 
A more detailed description of the apparatus can be found in Ref. [B7 . 



4 Data analysis 



In the offline analysis, cuts in the corrected coincidence time spectrum were 
applied to identify real coincidences and to eliminate the background of ac- 
cidental coincidences. A cut in the energy deposited in the first scintillator 
plane in the pion spectrometer was used to discriminate charged pions against 
protons. The Cerenkov signal was used to identify electrons in the electron 
spectrometer and to veto against positrons in the pion spectrometer. 

The true coincidences were observed in the peak of the accumulated missing 
mass distribution (E^ iss — IPmisgl 2 ) 1 ' 2 — M n , using an event-by-event recon- 
struction of (£ miss ,p miss ) = (u + M p - E n , q - pj. 

The cross sections were subsequently corrected for detector and coincidence 
inefficiencies (between +1.7% and +3.2%) and dead-time losses (between 
+1.4% and +4.6%). The detector efficiencies and their uncertainties were 
measured by the three-detector method, while the coincidence efficiency of 
the setup and the uncertainty of the dead-time measurement were determined 
in simultaneous single-arm and coincidence measurements of elastic p(e, e'p), 
which was also used to check on the acceptance for the extended targets. 



6 



The accepted phase space was determined by a Monte Carlo simulation which 
provided the (event-wise) Lorentz transformation to the CM system and in- 
corporated radiative corrections and ionisation losses of the incoming and 
scattered electrons and pions. Full track was kept of the particles' trajectories 
and their lengths in target and detector materials. Since exact energy losses 
are not known event- wise, we used the most probable energy losses for the 
subsequent energy loss correction in both simulation and analysis programs, 
and compared the corresponding missing mass spectra. The uncertainty esti- 
mates were based on relative variations of their content in dependence of the 
cut-off energy along the radiative tail. 

The uncertainty of the integrated luminosity originates only in the target den- 
sity changes due to temperature fluctuations within the target cell, while the 
electron beam current is virtually exactly known. Finally, a computer simula- 
tion was used to determine the correction factors due to the pion decaying in 
flight from the interaction point to the scintillation detectors, taking into ac- 
count the muon contamination at the target (correction factors ranging from 
x2.23 to x2.89 in different settings). The systematical errors of the pion de- 
cay correction factors were estimated from the statistical fluctuations of the 
back-traced muon contamination at the target. 



5 Results and discussion 



In the Born approximation, the coincidence cross section for pion electropro- 
duction can be factorised as 



23,28 



da 



dEL dQL dttt 



, dcr v 

v dm 



(5) 



where T v is the virtual photon flux and da v / dQ* is the virtual photon cross 
section in the CM frame of the final 7rN system. It can be further decomposed 
into transverse, longitudinal and two interference parts, 

+ £ L-T7^ + ^et(l + e)-— cos^ + e— — cos 2^ (6) 



dm dm ' u dm v ^ dm " d^: 



with the transverse (e) and longitudinal (e£ = Q 2 e/u>* 2 ) polarisations of the 
virtual photon fixed by the electron kinematics. The measured cross sections 
are listed in Table 0. 



In parallel kinematics (8* = 9^ = 0°) the interference parts vanish due to 
their sin#* and sin 2 ^* dependence. At constant Q 2 , the transverse and the 



7 



Table 2 

Measured cross sections in the p(e, e'7r + )n reaction. 



Setting (e) da/ dQ* Stat, error Syst. error 
[//b/sr] [jub/sr] [/xb/sr] 






834 


11 


14 





08 


(0.7%) 





41 


(3 


7%) 





500 


8 


40 





11 


(1.3%) 





31 


(3 


7%) 





219 


5 


96 





14 


(2.3%) 





19 


(3 


2%) 





742 


7 


48 





11 


(1.5%) 





18 


(2 


4%) 





437 


5 


36 





09 


(1.7%) 





10 


(1 


8%) 





229 


4 


55 





10 


(2.1%) 





07 


(1 


6%) 





648 


5 


03 





05 


(1.1%) 





08 


(1 


6%) 





457 


4 


19 





06 


(1.3%) 





08 


(1 


8%) 





259 


3 


46 





05 


(1.4%) 





07 


(1 


9%) 



longitudinal cross sections can therefore be separated using the Rosenbluth 
method by varying e, 

da v _ da T Q^_ da L 

dni ~ dn* lo* 2 dn*' u 

Not only the statistical, but also the e-correlated systematical uncertainties 
of the data were considered in our fit, whereas the e- independent systemati- 
cal errors were included in the final uncertainty of dor and da-^. The results 
are shown in Fig. || and Table ||. The Q 2 -dependence of the separated trans- 
verse and longitudinal cross sections can be seen in Fig. f|, together with the 
theoretical fits used to extract the form factors. 



Since values of v = —Q 2 /M 2 and W in this experiment were too high for 
a direct application of xPT, an effective Lagrangian model p3| , p4]| was used 
to analyse the measured Q 2 -dependence of the cross section, and to extract 
the nucleon axial and pion charge form factors. In the energy region of our 
experiment, the pseudovector 7rNN coupling evaluated at tree- level provided 
an adequate description of the reaction cross section. We included the s- 
and w-channel nucleon pole terms containing electric and magnetic Sachs 
nucleon form factors of the well-known dipole form with a 'cut-off' = 
0.843 GeV, the t-channel pion pole term with a monopole form factor F n (Q 2 ) = 
1/(1 + Q 2 /Al), the contact (seagull) term with the axial dipole form factor 
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Fig. 3. Least-squares straight-line fits to the measured cross sections for p(e,e'7r + )n 
at W = 1125 MeV, for three values of Q 2 . The smaller error bars correspond to 
statistical, the larger ones to the sum of statistical and systematical errors. 

Table 3 

The results of the L/T separation based on the least-squares fit to the data. 



Setting (Q 2 ) 
[GeV 2 /c 2 ] 



da T /dn* 

[/ib/sr] 



(Q 2 /^ 2 ) daJdQ* 
[/ib/sr] 



0.117 4.160 ±0.165 sta t±0.202 sys 8.394 ± 0.254 stat ± 0.48 l sys 
0.195 3.080 ± 0.139stat ± 0.051 sys 5.747 ± 0.284 stat ± 0.246 sys 
0.273 2.406 ± 0.088 stat ± 0.056 sys 4.012 ± 0.187 stat ± 0.052 sys 



Ga(Q 2 ) = Ga(0)/(1 + Q 2 /MD 2 , and the s-channel A-resonance term. Vector 
meson exchange contributions in the t-channel were found to play a negligible 
role in the charged pion channel. 

Due to cancellations between higher partial waves and interference terms with 
the s-wave, dor is predominantly sensitive to the E 0+ (mv + ) amplitude and 
therefore to Ma. On the other hand, the pion charge form factor appears 
in the longitudinal amplitude L +(n7r + ) only at order (9(/i 2 ,/«/), and the s- 
wave contribution to dcr L amounts to 10 % only. Due to the contributions of 
the higher partial waves, however, the longitudinal cross section di7 L is quite 
sensitive to A^, which is a bonus 'by-product' of the analysis. 

Since current conservation is violated if arbitrary form factors are included, 
gauge invariance was imposed by additional gauge terms in the hadronic cur- 
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Fig. 4. Separated transverse and longitudinal cross sections. The solid line shows 
our fit with da T /dQ*(Q 2 = 0) fixed to (7.22 ± 0.36) ^b/sr; the dotted line is the 
unconstrained fit. In the longitudinal part the fits are almost indistinguishable. The 
smaller error bars correspond to statistical, the larger ones to the sum of statistical 
and systematical errors. 

rent. These terms modify the longitudinal part of the cross section and there- 
fore influence the pion pole term and thus the extracted value of A w . However, 
this procedure does not affect the transverse part of the cross section. The 
axial mass can then be determined from the <3 2 -dependence of dor, and the 
result can be compared to the prediction of In the theoretical fit, the 

transverse part is fitted first. 
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We have used two different techniques to determine the axial mass: 



(I) Only the axial mass was varied, whereas the value of dor at Q 2 = 
was fixed by extrapolating the transverse cross section (i. e. the Eo + (mr + ) 
amplitude) to the value of the photoproduction angular distribution at 9 n = 
0°. A precise cross section at the photon point (Q 2 = 0) is very helpful for 
our analysis and can reduce the uncertainty in determining the axial mass 
considerably. Unfortunately, there are no experimental data available around 
W = 1125 MeV and forward angles. The only measurements in this energy 
region are performed at pion angles larger than 60° with large deviations 
among the different data sets. Therefore we have used a value at the photon 
point obtained from different partial-wave analyses of the VPI group p9|| and 



of the Mainz dispersion analysis ||30|| . This yields to a weighted-average cross 
section at the photon point of (7.22 ± 0.36) /ib/sr and this value was used as 
an additional data point. The corresponding value of Eo + (mr + ) is also well 



supported by the studies of the GDH sum rule ]31j] and by the low energy 
theorem (Kroll-Ruderman limit). 

(II) The three data points alone were fitted, while the value of the transverse 
cross section at Q 2 = was taken as an additional parameter, with the result 
da T (0) = (7.06 ± 1.12) /ib/sr. 



The best-fit parameters for the transverse cross section were then used to fit 
the longitudinal part. 

Using the first and preferred procedure, we find from the transverse cross sec- 
tion M A = (1.077±0.039) GeV, corresponding to (ri) 1 / 2 = (0.635±0.023) fm. 
From the longitudinal part, we obtain = (0.654 ± 0.027) GeV, correspond- 
ing to (r 2 ) 1 / 2 = (0.740 ± 0.031) fm. The second procedure leads to the follow- 
ing results: M A = (1.089 ± 0.106) GeV or {r\) l l 2 = (0.628 ± 0.061) fm, and 
A^ = (0.658 ± 0.028) GeV or (r 2 ) 1 / 2 = (0.734 ± 0.031) fm. 



6 Summary and conclusions 



We have measured the electroproduction of positive pions on protons at the 
invariant mass of W = 1125 MeV, and at four- momentum transfers of Q 2 = 
0.195 (GeV/c) 2 and 0.273 (GeV/c) 2 . In conjunction with our previous mea- 
surement at Q 2 = 0.117 (GeV/c) 2 |23H , we were then able to study the (in- 
dependence of the transverse and longitudinal cross sections, separated by the 
Rosenbluth technique for each Q 2 . The statistical uncertainties were between 
0.7% and 2.3 %, an improvement of an order of magnitude over the result of 
Ref. fl33"| . The systematical uncertainties were estimated to be between 1.6% 
and 3.7 %, and are expected to decrease significantly in the future experiments. 



11 



We have extracted the axial mass parameter Ma of the nucleoli axial form fac- 
tor from our pion electroproduction data using an effective Lagrangian model 
with pseudovector 7rNN coupling. Our extracted value of Ma = (1.077 ± 
0.039) GeV is (0.051 ± 0.044) GeV larger than the axial mass M A = (1.026 ± 
0.021) GeV known from neutrino scattering experiments. Our result essentially 
confirms with the scaled-error weighted average Ma = (1.068 ± 0.017) GeV of 
older pion electroproduction experiments. If we include our value into the 
database, the weighted average increases to Ma = (1.069 ± 0.016) GeV, and 
the 'axial mass discrepancy' becomes AMa = (0.043 ±0.026) GeV. This value 
of AMa is in agreement with the prediction derived from x?T, AMa — 
0.056 GeV. We conclude that the puzzle of seemingly different axial radii as 
extracted from pion electroproduction and neutrino scattering can be resolved 
by pion loop corrections to the former process, and that the size of the pre- 
dicted corrections is confirmed by our experiment. 

Theoretical input needed to extract the axial mass and the pion radius has 
naturally led to some model dependence of the results. The dominant contri- 
bution to pion electroproduction at W = 1125 MeV is due to the Born terms. 
These are based on very fundamental grounds and the couplings are very well 
known. For our purpose the electric and magnetic form factors of the nucleons 
are also accurately known whereas the remaining two (the axial and the pion 
form factors) are the subjects of our analysis. In parallel kinematics we are in 
the ideal situation where the pion form factor contributes only to the longitu- 
dinal cross section, therefore reducing very strongly the model dependence of 
the transverse cross section and consequently of the determination of the axial 
mass. Furthermore, the A resonance, which plays the second important role in 
our theoretical description, also contributes with a well-known Ml excitation 
to the transverse cross section, while the longitudinal C2 excitation gives rise 
to a larger uncertainty due to the less known C2 form factor, currently under 
investigation at different laboratories. 

Altogether, the model dependence is smaller for extracting the axial mass 
than for the pion radius. This could partly explain the discrepancy in the 
different values of the pion radius between our analysis and the analysis of pion 



scattering off atomic electrons [32] . However, as in the case of the axial mass, 
an additional correction of the pion radius obtained from electroproduction 
experiments is very likely. This should be investigated in future studies of 
chiral perturbation theory. 
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